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Abstract. A simple proof of Ramanujan's formula for the Fourier transform 
of |r(a + is given where a is fixed and has positive real part and t is real. 
The result is extended to other values of a by solving an inhomogeneous ODE 
and we use it to calculate the jump across the imaginary axis. 



1. Introduction 

The behavior of the modulus of the gamma function along vertical lines plays 
an important role in analytic number theory pp 767-775].) Whenever the real 
number a is not a negative integer or zero, the function 

t ^ \r{a + it)f 

is an even function in the Schwartz class, i.e., it decays rapidly at infinity along 
with all derivatives. When a > 0, one can prove the following estimate due to Lerch 
m P- 15]): 



Va^ -I- V sinhTTt' 

where A satisfies 1 < A < vT+i^, and if a = ^ , we have the closed form expression 

2 



TT 



cosh nt 



It is natural to ask what we can say about the Fourier transform of the Schwartz 
class function t f->- |r(a -f it)|^. The following beautiful formula was proved by 
Srinivasa Ramanujan in jlO^. For a > 0, we have 

(1.1) 



J \r{a + it)fexp{-i^t)dt^VTTr{a)r(^a+^^ ^cosh 
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This gives the Fourier transform of the rapidly decreasing (i.e., Schwartz class) 
function t i— > |r(a + it)| , and by Fourier Inversion, also the Fourier transform of 

^ I— >■ ^cosh|^ . The existence of either of these transforms as a closed form 
expression is a minor miracle. Setting a = ^ in (II. ip gives the well-known formula 

, cosh nt cosh | 

revealing the close relation with the Gaussian function. 

The integral (11.11) is the simplest in a class of integrals considered by Mellin, 
Ramanujan, Hardy, Hecke and others (0] p. 98, pp. 103-104, and p. 492]). Mellin 
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studied integrals of this form in connection with representations of solutions of first 
order difference equations ([71 pp. 344-345]) and solutions of (generalized) hyper- 
geometric equations pp 148-149]). Integrals of this form also play an important 
role in analytic number theory ([5j p. 152]) and in the theory of real quadratic 
number fields ([9, p. 58]). 

In this note we consider what happens to the formula ()1.1|) when — 1 < Re (a) < 0. 
It is clear that the formula as it stands must fail when Re(a) is negative and a ^ Z 
since for such a, the function 1 1— > ]r(a-|-ii)]^ is rapidly decreasing, but ^ cosh^ ° ^ 
on the right hand side blows up exponentially. We obtain, using Fourier analysis, an 
inhomogeneous ordinary differential equation for the integral when — 1 < Re (a) < 0. 
The formula assumes a very different form when Re(a) < 0. To make the paper 
self-contained we offer using techniques of Fourier analysis, an alternate proof of 
(frTj) when Re (a) > 0. 

We will prove Ramanujan's formula in the following form, where we allow a to 
be complex: 

Theorem 1.1. When Re(a) > 0, 

27rr(2a) 



/OO 
r(a -I- it)r{a - it) exp{-i£,t)dt 
-CO 



4*^ cosh^' 



°(« 



When a is real, the left hand side of (II. 2p reduces to that of (|1.1[) . since r(z) = 
r(z) for all z in the domain of f, and therefore 

r(a + it)r{a - it) = \r{a + it)\'^ , 

and in the right hand side, the two expressions coincide, thanks to the classical 
duplication formula for the gamma function ([HI p. 35], [131 formula 3.8]): 

(1.3) v^r(2a) = 22"-ir(a)r (^a+ i 

If we try to prove Theorem 1 1.1 1 bv formally computing the Fourier transform, we 
soon run into a familiar obstacle, namely, the appearance of an oscillatory integral 

exp(i^(a; — y))dx. 

The standard way to cope with this (see [TH p. 37]) is to introduce a Gaussian 
convergence factor exp(— ea;^) (with e > 0) in the integrand, compute the integral, 
and let let e ^ 0+. This is done in Section l2] below. 

We then ask how to modify (|1.2p when Re(a) is negative. We show that in 
Section |31 the question can be reduced to the solution of an inhomogeneous linear 
ODE, whose solution can be reduced to (repeated) quadratures. When —1 < 
Re(a) < we prove the following: 

Theorem 1.2. When Re(a) G (—1,0), we have 

r(a + it)r{a — it) exp{—i£^t)dt = 



, , ^, , e'^'' , /•« cosh(a(e-s))^ ' 

coshlaf ) / „ . „ , ds + / „ , „ , ^ ds 

^ ^' 'o cosh"'+^ (f) Jo cosh2'^+2 (f ) 
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As a consequence, we find the jump in the Fourier transform of 1 1— >■ r{a—it)r{a+ 
it) as a crosses the imaginary axis in Corollary 14. II 

Theorem 11.21 shows that when a < 0, the Fourier transform is no longer in closed 
form, but can still be expressed in terms of quadratures. It will be seen in the proof 
that this method can be extended to obtain similar formulas for vertical strips of 
the form — (n + 1) < Re(a) < — n, for any positive integer n in terms of multiple 
quadratures. 

2. Proof of Theorem 11.11 

We begin with a useful formula that goes back to Binet ('2', p. 136]). We recall 
that the beta function B{p,q) is defined for Re{p) > and Re(g) > to be 



B(p,(?)= f\''"\l-s)''-'ds, 
Jo 



and the basic relation between the beta and gamma functions ([TTl p. 35], [131 p. 
313]): 

(2.1) r{p)r{q)^r{p + q)B{p,q). 

Lemma 2.1 (Binet). When Re{p) > and Re{q) > 0, 

r°° p{p-q)u I p{q-p)u 
(2-2) B(|,,g)=/ —-^----d.. 



(e" + e-")P+9 

Proof. We set s ~ (1 + e^")~^ in the integral defining B(p, q): 



/oo 



g2tip _|_ ^2uq 



du 



du, 



, (1 + e2")P+9 

since the beta function is symmetric in p and q. Multiplying the numerator and 
denominator of the integrand in the last integral by exp {—u{p + q)), we get (|2.2I) . 

□ 



We now prove Theorem 11.11 By virtue of the beta-gamma relation (12.1 
r(a + it)V{a - it) = B(a + it,a- it)\'{2a). 
Invoking Lemma [531 we see that the left hand side 

/oo 
r(a + it)r{a — it) cxp{—i£_t)dt 
-oo 

of (|1.2p is equal to 

(e" + e"") 



r(2a) / exp(-^^0 ( / W^T^^^" ) dt 



dt. 



POO f / i^oo / 2itu , -2itu\ \^ 

(2.3) = r(2a) ^hm |exp(-et^ - .^0 \. ^ ^^.^J du) | 

Where in the last line, we have inserted a convergence factor exp(— et^) (with e > 0). 
Note that the function 

/■oo g2itu g — 2itu 

t^ I -. rrr- du 

, (e" + e-")2'' 
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is rapidly decreasing, since it can be interpreted as the sum of Fourier transforms 
of rapidly decreasing functions, and therefore is in particular in L^. Therefore, we 
are justified in applying the dominated convergence theorem, and writing (|2.3p as 



r(2a) lim 



where in the last line the use of Fubini's theorem to interchange the order of integra- 
tion is easily justified, since 1 1— > exp(— ei^) is integrable. Inserting the well-known 
formula ([HI p. 41]) 



/■°° / 2 ^ , ^ c2 

/ expl— et — ictjdt = \ — exp 

J-oo ye V 4e 



for the Fourier transform of the Gaussian into each term of the second integral 

factor, we get 

(2.4) 

r /-oo e-M^ 1 e"^^ 1 1 

-v/7rr(2a) lim < / — • z^du -f / — • z^rdu > 

'^-^^{j-oo (e" + e-")^" 7-00 (e"-He-«)2'» J 

The change of variables w = | + yfe w transforms the first integral in p.4p into 

(2.5) / J . 

J-oo (eHv^™ + e-5-V^«')2a 

Since the denominator is always greater than or equal to 1, in this integral, the 
integrand is dominated by the L function w i— >■ e"™ . We may therefore apply the 
dominated convergence theorem and conclude that the integral (|2.5p converges as 
e 0+ to 

r- 1 



(2.6) 



[el+e-^y 4'^cosh-(|) 



The change of variables v — —u shows that the second integral in (|2.4p is equal 
to the first for each e > 0, and therefore, its limit as e — 0+ is also equal (|2.6I) . 
and therefore, the upshot of our computation is that the expression in (j2.4p . and 
therefore, the Fourier transform on the left hand side of (II. 2p is equal to 

27rr(2a) 



4° cosh^" (l) 



which is the required right hand side. 

Remark: That the first integral in (|2.4p converges as e — > 0+ to the quantity 
in (|2.6p may also be deduced from the following fact: for a fixed ^, as e — > 0+, the 

1 — (2u — 

first factor u i— > e^^e ~^ of the first integral in (|2.4I) converges in the sense of 
distributions to \/TTdi, where for p G R, we denote by 6p the delta distribution on 
M supported at the point p. 
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3. Proof of Theorem 11.21 
Suppose that Re(a) is not a negative integer or 0. For real ^, we define 

/oo 
ar{a-it)r{a + it)e-'*^dt. 
-oo 

This is well-defined, since the integrand is easily seen to be a rapidly decreasing 
function of t. When Re(a) > 0, Theorem 11.11 along with the duphcation formula 
(|1.3p shows that 

(3.2) /(a, = a«'(a) sech^" 
where we denote 

and sech(<) = (cosh(t)) ^ . The basic property of the function I{a, ^) is given in the 
following lemma: 

Lemma 3.1. //Re(a) is not a negative integer or zero, we have 

(3.3) ($-'^0'^"'^^=(^)'^" + ''^^- 

Proof. Using the functional relation r(z) ~ — — — — - twice, the right hand side of 

z 

p.ip can be rewritten as 

(3.4) I{a,^)^J (^-^^^r{l + a~it)V{l + a + it)e-'*^dt. 

To get rid of the factor + in the denominator in the integrand in (13. 4p we 
differentiate under the integral twice with respect to ^, and subtract a^I{a,^). We 
obtain 



di 



2 



I{a, = y" r(l + a - it)r{l + a + it)e-''^dt 

i ^^p ] r(l + a - it)r{l + a + it)e~'*^dt 

/OO 
r(l + a - it)r{l + a + it)e~'^^dt 
-OO 

(^TTt) / (« + l)'"(l + a-«Or(l + a + j<)e""«dt 
/(a + 1,0. 



—a 
a + 1 



□ 



We now consider the situation in which Re(a) > —1, and Re(a) ^ 0. Then 
Re(a + 1) > 0, and consequently, putting the value of I{a + l,^) from (13. 2p into 
(13. 3|) we obtain that 

(;| - = (^) + 1) (f) . 
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We have thus obtained the foUowing inhomogeneous linear differential equation 
satisfied by the function ^ !-> I{a,£_): 



(3.5) - a'j /(a, ?) = -a^a + l)sech'^+'^. 
The general solution of the associated homogeneous equation 

is ^ I— Ae'^^ + Be'""^ , where A and B are constants. We proceed to find a particular 
solution of p.Sp by the standard method of variation of parameters. We seek the 
particular integral in the form 

(3.6) e^^^i(C)e'^«+«2(e)e-"«, 

where the functions vi and V2 are determined by the system of equations 

«i(^)e°^ - «2(0e"''* = - $(a + l)sech2 ^ 



2 



Solving the above pair of equations for v'l and , representing vi and f 2 as indefinite 
integrals, and substituting in (13.61) gives the particular integral 

The complete solution of (13. Sp is thus 

The constants A and i? are to be determined by side conditions such as initial or 
boundary conditions. The Riemann Lebesgue lemma (14, p. 106 fi']) which asserts 
that the Fourier transform of an integrable function decays at infinity, provides the 
requisite boundary conditions 

(3.8) lim/(a,O=0, lim l{a,^)^0 

Taking the limits as ^ — >■ 00 and ^ — >■ — 00 in p.7p . and using p.Sp we get the 
following values of A and B, we see that if Re(a) < 0, we have 



(3.7) 



(3.9) A = B = -^^ rsech-+^(0e- 



'ds. 
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For use in the next section, we note that when Re(a) > 0, the same method gives 
us that 

(3.10) A^B^ r sech-+^ (^) e^ds. 

Substitution of the values of A and B obtained in (|3.9I) in (|3.7p . and simphfying 
the expression, gives us that for —1 < Re(a) < 0, we have 

/(a, ^) = _$(a + 1) I jJI^ scch2"+2 (0 cosh (a(^ - s)) ds+ 

cosha^ scch2°+2 (^0 e^^'ds 

and this is equivalent to Theorem 1 1.2 1 in view of the definitions of I{a, ^) and $(a). 

Remark: Note that the method of proof can be used to generate formulas for 
/(a, £_) (and consequently the left hand side of (|1.4p ') for Re(a) < — 1 by a recursive 
process. These formulas will involve repeated quadratures. For example, suppose 
for a positive integer n that /(a, ^) has been determined for Re(a) S (— —n + 1). 
Then the right hand of (13.31) is known when Re(a) G (— n — l,n). Consequently, 
we have a linear inhomogeneous ordinary differential equation to find /(a, ^) when 
Re(a) G (— n — 1, n), and this can again be solved using the method of variation of 
parameters, yielding a solution in terms of multiple quadratures, where the inte- 
grands now involve only elementary functions. 



4. The Jump Across the Imaginary Axis 

As we have seen in the previous section, for fixed real ^, the left hand side of (|1.2p 
is a holomorphic function of a, provided Re(a) is not a negative integer or zero, and 
the technique of the last section may be used to obtain formulas representing this 
function in terms of integrals of elementary functions. Therefore we can find the 
jump of this function of a along the lines {Re(a) = — n}, where n is a non- negative 
integer. We now carry out this program for the simplest case of n = 0. Denote for 
convenience 

/oo 
r{a-it)r{a + it)e-'*^dt, 
-OO 

and for a point p on the imaginary axis consider the jump 

RG(a)>0 Ro(a)<0 

Corollary 4.1. For each real ^, and for purely imaginary p ^ 0, we have 

[[J(p,e)]] =4^cosh(pe)r(2p). 
Proof. Notice that, the function J above and the function / of (|3.ip are related by 

Set p = ic, where c G R, and let [[/(ic, ^)]] denote the jump across the imaginary 
axis: 

[[I{ic,m= lini l{a,0- lim /(a,e). 

a — >-ic a—^ic 
Rc(a)>0 Re(a)<0 
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In order to prove the result, it suffices to sfiow tliat for every nonzero real c, we 
have 

[[I(ic, ^)]] = 27r cos cCr(l + 2ic). 
The integrals in p.7p are continuous functions of a since the real part of a + 1 is 
positive. Hence contribution to the jump of /(•, ^) at the point ic on the imaginary 
axis arises from the discontinuity in the coefficients A and B. From equations (I3.9P 
and (|3.10p . the difference in the right and left hand limits of Ae°-^ is 

Using Binet's formula with p — 2ic+ 1 and q ^ I wc simplify the last expression to 

2i+2*<=$(l + ic)e*'=«B(l + 2ic, 1) = 2^"'V^r{l + ic)r{ic + i)e''=«. 

Using the duplication formula (|1.3p to simplify the last expression, we see that the 
difference in the right and left hand limits of Ae"-^ at ic is 

7re*"«r(l + 2ic). 

Similarly, the difference in the right and left hand limits of Be^""^ at ic is 

7re-'"«r(l + 2ic). 

Thus we have 

[[/(ic,e)]] = TTie'"^ + e-*"«)r(l + 2ic) = 27rcosc^r(l + 2ic) 
as required. □ 
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